Exact results for the Spectra of Bosons and Fermions 

with Contact Interaction 

Stefan Mashkevich^ 

Schrodinger, 120 West 45th St., New York, NY 10036, USA 

Sergey Matveenko^ 

Landau Institute for Theoretical Physics, Kosygina Str. 2, 119334, Moscow, Russia 

Stephane Ouvry^ 

Laboratoire de Physique Theorique et Modeles Statistiques"^ 
Bat. 100, Universite Paris-Sud, 91405 Orsay, France 

Abstract 

An A^-body bosonic model with delta-contact interactions projected on the lowest 
Landau level is considered. For a given number of particles in a given angular momentum 
sector, any energy level can be obtained exactly by means of diagonalizing a finite matrix: 
they are roots of algebraic equations. A complete solution of the three-body problem 
is presented, some general properties of the A^-body spectrum are pointed out, and a 
number of novel exact analytic eigenstates are obtained. The FQHE A^-fermion model 
with Laplacian-delta interactions is also considered along the same lines of analysis. New 
exact eigenstates are proposed, along with the Slater determinant, whose eigenvalues are 
shown to be related to Catalan numbers. 

PACS numbers: 71.10-w, 71.70Di, 05.30.Jp 



1 Introduction 

There has recently been considerable interest in Bose-Einstein condensates, following their 
experimental discovery in atomic vapors |21 E] ■ Particular attention has been devoted 
to fast-rotating condensates, since such systems were demonstrated experimentally to 
form vortices jH |H] , like in superfluid ^He; systems with ever bigger numbers of vortices 
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were subsequently observed [3 13 E] • A plausible theoretical model describing the system 
is that of two-dimensional (perpendicular to the axis of rotation) trapped bosons in a 
harmonic well with delta-function repulsion. The fast rotation implies that the particles 
are in a strong effective magnetic field, thereby allowing for a lowest Landau level (LLL) 
analysis. States with a given number of vortices correspond to the sector of fixed angular 
momentum L in the A^-body problem. 

Apart from its application to Bose-Einstein condensates, the problem is interesting 
per se, as an A^-body quantum-mechanical model which can be solved exactly, at least in 
part. This model has a clear relationship with another well-known two-dimensional setup 
in which the classical interaction force is also zero at nonzero distance — the anyon model 
p. In the latter, only a small subset of the spectrum can be found exactly [TU], whereas 
within the lowest Landau level, all eigenstates and eigenvalues are explicitly known [TTj . 
One might expect that the delta-interaction LLL problem at hand would be solvable in 
the same way, too. But this does not appear to be the case. An obvious explanation for 
this difficulty lies in the very definition of the model itself: delta interactions are known 
to be ill-defined in 2 dimensions, since the delta term is too singular with respect to the 
kinetic term. This singular interaction has to be regularized in some way. In fact, the 
anyon model itself can be viewed as a regularized version of the delta-interaction model: 
the long-distance quantum Aharonov-Bohm interaction, at the heart of the anyonic in- 
termediate statistics, smoothes out the perturbative divergences due to the short-range 
contact interactions also present in the anyon model (exclusion of the diagonal of the 
configuration space). 

Here, in the model considered, the regularization used is the LLL projection itself. It 
can be checked that restricting to the LLL subspace allows for a finite contact interaction 
to all orders in perturbation theory. 

It was conjectured, based on numerical results, that the ground state for L > N has a 
simple analytical form [T2]; it was proved later [Tnilllj that the state in question is indeed 
an eigenstate. A hierarchy of A^-body exact eigenstates for L = 4 was discovered J3]. The 
other results are mainly numerical: the collective vs single-particle nature of the excited 
states was explored, good single-particle approximations in the L <^ N regime discovered 
^14j : very recently, a composite- fermion ansatz was suggested, and it was shown that the 
CF wave functions at L = A^ become exact in the thermodynamic limit |15j . 

In this paper, we concentrate on exact results. In Sec. 2, we formulate the problem, 
perform the projection onto the lowest Landau level, and show that the interaction can 
be diagonalized separately in each given angular momentum sector. This means that any 
level can be obtained by diagonalizing a finite matrix; no basis truncation is necessary. In 
Sec. 3, we discuss some general properties of the spectrum and review a number of cases 
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where the eigenstates have a simple polynomial form and the corresponding eigenvalues 
are rational numbers. In the next section it is shown, by direct counting, that those 
simple eigenvalues exhaust the whole spectrum of the three-body problem, which is thus 
completely solved. Section 5 deals with the A^-body problem. Given the nature of the 
wave functions involved, it is convenient to consider states with a given value of L and 
varying A^. At L = 2, 3, all the A^-body eigenstates fall into the category of the simple 
states of Sec. 3. For each of L = 4 and L = 5, apart from the simple ones, there is one 
extra A^-body eigenstate: for L = 4, it is the one found in Ref. Jl]; for L = 5, it is new. 
Starting with L = 6, most eigenvalues are irrational, but all of them are roots of algebraic 
equations. An algorithm is built that lets one find them all. 

Finally, in Sec. 6, a fermionic model introduced so that the Laughlin FQHE wave- 
functions are its ground states is considered, again in the LLL. Via a mapping onto 
an equivalent bosonic model, it can be analyzed following the same approach as in the 
bosonic case. A few eigenstates are given as illustrations of the procedure. The A^-body 
Slater determinants, with any A^, have eigenvalues which are directly related to Catalan 
numbers. 



2 Formulation and method of solution 

Consider A^ bosons in two dimensions in a harmonic well with pairwise contact interaction: 
V = VJ2i<j (^(^i — Tj). In complex coordinates {z = x + iy), the Hamiltonian is 

N _ 2 N N 

H = ~2j2dA + — J2z,z, + V S'iz,-z,). (1) 

i=l i=l i<j = ^ 

Projection on the lowest Landau level. As said above, the delta interaction is too 
singular in two dimensions and as such is not properly defined. One way to regularize 
the Hamiltonian is by projecting it on the "lowest Landau level" (LLL) subspace. This 
means that upon extracting the long-distance harmonic damping exponential factor one 
restricts to the set of A^-body wave functions 

^(ri, ra, . . . , Vn) = exp ^-^ ^ ZiZ^^ f{zi, Z2, Zn) , (2) 

where f{zi,Z2, .--^Z]^) is analytic. 

Indeed, if a magnetic field were added, the single-particle eigenstates corresponding 
to the LLL (Landau level number n = 0, angular momentum / > 0) would be in the 
LLL-harmonic eigenstates basis 
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{z,z\Q,l)=l^-^j ^'e-^--", (3) 

where Ut = ^Juf+^, Uc being half the cyclotron frequency. Coming back to the 
purely harmonic problem, consider now the projector on the "LLL" Hilbert space Pq — 

E£o|o,0(o,/| 

{z, z\Po\z\ z!) = ^^e-t(-+-'-"'-2-') (4) 

71 

A generic state belonging to the LLL, = J2i1o'^i\^j^)j analytic up to the Gaussian 
factor, 

{z\^/j) = f{z)e-^^' (5) 

with 

oo 

/(^)=E«'^'- (6) 

1=0 

Projecting a one-body Hamiltonian on the LLL amounts to 

{z, z\PoHPo\4^) = {z, z\PoH\il;) = / dz'dz'-e-^^''+''''-^'''^H{z')f{z')e-^'''' , (7) 

J n 

that is to say, for a one-body potential V{z, z), to the eigenvalue equation 

/ dz'dz'-e-'^^''''-'''^ [a; + ujz'd' + V{z', z')] f{z') = Ef{z) . (8) 

In the A'"-body case with V — VJ2i<j ^"^{^i — ^j) one obtains a product of N integrals 
of /(^i, zjv). Using the Bargmann identity 

^/ci.'ci.-'e-(^'^-'-')M.') = M-), (9) 

TT J 

after integrations, the projected Hamiltonian reads (for simplicity, a; = 1 and V/ (27r) = 1) 

Hj^LL^L + V, (10) 
where the angular momentum operator is 

AT 

L^Y^^i^i (11) 
1=1 

and the LLL-projected contact interaction 

N 

y-Y. P^3 , (12) 
i<j=l 
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Pij being the operator^ that replaces both Zi and Zj in the wave function with the coor- 
dinate of the center of mass of the pair (ij), keeping all the other coordinates intact: 

P.,/(.. .,.„...,.,,...)=/(.. .,^,...,^,...) . (13) 

A generic bosonic eigenstate of L is a symmetrized monomial in the Zj's: 

fh...lN = ^ 

with 

N 



■ N 
i=l 



(14) 



L = Y.k, li>k+i>0. (15) 



i=l 

Here the symmetrization operator acts as 

S[f{zi,z2,. . . ,zn)] = ^^/(2;7'i,2;p2,...,2p^) , (16) 

the sum being over all A^! permutations P of (1, ... , A^). 

The multiplicity G{N, L) of a level with given values of A^ and L is the number of 
ways that L can be represented as a sum (fT^ . i.e., the number of unordered partitions 
of L into A^ nonnegative addends. [E.g., for A^ = 3, L = 6, the (7(3, 6) = 7 partitions are 
(6,0,0), (5,1,0), (4,2,0), (4,1,1), (3,3,0), (3,2,1), (2,2,2).] It is directly connected to 
a well-studied quantity p]\f{M), the number of unordered partitions of M into A^ positive 
addends: By adding 1 to each addend in each partition contributing to G{N, L), one sees 
that G{N,L) = pn{L + N). This is a polynomial of degree (A^ - 1) in L: G{N,L) = 
L^-^/[N\{N - 1)!] + 0{L^~^). The generating function is 

oo -| oo 

n— -^ = 1+ E G{N,M-N)x'^^z^ . (17) 

k=l ^ ^ M,N=1 



^For V = V J2i<j<k=i ^'^i^i ^ ^ ^k), when the Pfafhan states are known to play a central role, 

one would have in the LLL 

N 

i<j<fc=i 

where Pijk is the operator that replaces Zi, Zj and Zk in the wave function f (. . . , Zi, . . . , Zj, . . . , Zk, ■ ■ .) 
with the coordinate of the center of mass of the triplet (ijk), keeping all the other coordinates intact: 

p r( \-ff ^'^ + Zj + Zfc Zj + Zj + Zfc \ 

-t^ijkj [■ ■ ■ , Zi, . . . , Zj , . . . , Zk, ■ ■ ■) — / I • ■ • J „ J ■ • ■ J Q I ■ • ■ J Q J • ■ • I ■ 
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The problem reduces to the diagonahzation of V in the G{N, L)-dimensional subspace 
spanned by {fij^...iff}, for any given and L. A generic state befonging to that subspace 
is 

f{zi, Z2, Zn) = X! (^h-lNfh-lN (18) 
{h} 

with /ji.../^ defined by Eq. (dH), ai^...it^ arbitrary coefficients, and the sum being over all 
G{N,L) possible sets {Ik} defined by Eq. (fT3j) . The eigenvalue equation is 

Vf = Ef (19) 

[the total energy, according to (fTUj). is L + E]. One has 

N N 

vf= y: p^,f = E E P^jfh...i. , (20) 

and upon expanding all the powers of ^^-^ in Pijfi^...i,^, this becomes a polynomial in the 
Zi's. In Eq. (fTHj). with account for (IHj), ai-^...i^ is the coefficient at Hk^k- Hence, in the 
polynomial pn|) . the coefficient at Ha; -^if must be equated to Eai-^,,,i^, for each {/fc}; one 
obtains G{N,L) equations for as many coefficients a^^.../^. Every level can thus be found 
exactly by diagonahzation of a finite matrix, without the usual errors induced by basis 
truncation. 



3 General properties 

Center of mass. Introduce the center-of-mass coordinate 

1 ^ 



N . , 

1=1 



and the relative coordinates 



obviously. 



Si = Zi-Z ; (22) 



N 



E = • (23) 



i=l 



For any function of the center-of-mass coordinate, PijE{Z) = F{Z), hence 

VF{Z) = ^^^^'^^ F{Z) . (24) 
Moreover, for any eigenfunction f{zi, Z2, z^) of the interaction Hamiltonian, 

Vf{zi, Z2, zn) = Ef{zi, Z2, zn) , (25) 



also 

VFiZ)f{z,, Z2, zn) = EF{Z)f{z,, ^2, zm) (26) 

— a center-of-mass excitation does not change the interaction energy (although it does 
affect the total energy because the angular momentum changes). Thus, for each {N,L) 
state, there is a "tower" of CM excitations above it: {N,L + 1), {N,L + 2), etc. It is 
enough to find the "pure relative" states, whose number is 



G{N, L) = G(N, L) - G{N, L-l] 



(27) 



Known exact eigenstates. In the A^-body case, denote 



fh. 



S 



■ N 

.1=1 



{21 



One has 



N 



vh.., = E 

i<j=i 



Zi ~t~ Zj 



2 N 



N 



[Zi + z^ 



z \2 N 



N 



k=l 

Zi Zj 



4 i<j=i ZiZj ^^-^ 



N 



i<j=i 



/ k=l 



N 



N{N-1)-N+ J2 ^ 



N 



N{N-1)+ J2 ^ 



N 



k=l 



k=l 



[at the very last step, Eq. was used]. 
More generally, for 1 < L < A^, 



L 



f (iV-^-l)A....o...o 



(29) 



(30) 



(see Refs. \13\ IT^ for two different proofs). It has been conjectured numerically that this 
is the ground state for any L < N [12 . Obviously, for L = it reduces to (jSHj)- Note 
also that /10...0 = because of Eq. (|^. 

Further, in the 3-body case, for any L > 1, one has 



LOO 



zL 



i<j=l 
1 ^ 

- E 

i<j=l 



k=l 



1 + 



2 

2L-1 



Zi -\- Zj 

2~~ 

/loo 



z^ -z^ 



+ 3/. 



LOO 



(31) 
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since Zi + Zj = —Zk for i j k. 
In the 4-body case, for any odd L, 

^/iooo = 3/looo ; (32) 

since in this case Zi + zj = —z^. — zi for i ^ j ^ k ^ I. 

In the A^-body case, for L = 2 and L = 3, respectively, /20...0 and /30...0 are eigen- 

states of the form (jHn|) since, trivially, Y^k=i ^fc = ~ '^ki ^kZi and, less trivially, J2k=i ^1 = 
TV ~ — 

^k,l,m ^kZlZm- 

"Slater excitations". Consider the square of the A^-particle Slater determinant, 

N N 

s'n= Yi (^. - = n (^^ - • (33) 

i<j=l i<j=i 

Obviously, PijSl^f = for any f{zi, . . . , zn)- Hence, a "Slater excitation" of any function 
is an eigenfunction of V with zero eigenvalue: 

VS],^f{zi,Z2,...,ZN) = 0. (34) 

(The presence of a Slater determinant nullifies the probability for the positions of any 
two particles to coincide, hence the delta interaction has no effect.) An even power of 
the determinant is required in order to preserve the symmetry of the function. This is 
somewhat reminiscent of the idea of "composite fermions" J3] . 



4 Three-body problem 



For = 3, it turns out that the simple eigenstates described in the previous section 
exhaust the whole spectrum. 

The number of pure relative eigenstates for any L > is 



G(3,L) 



6 



+ 1 — (^L mod 6,1 5 (35) 



i.e., (5(3, 0) = 1, G(3, 1) = 0, G(3, 2) = . . . = G(3, 5) = 1, G(3, 6) = 2, and so on with 
period 6. Here is the systematics of these states. For any L 7^ 1, there is a state fioo with 
the eigenvalue 1 + (— 1)^^/2^"^, Eq. (jHT|) . These make up all the states for L < 5. [Recall 
that /no, Eq. (j^ . is the same as /2oo-] Now, upon any (3, L) state / there is a tower of 
(3, L + 6n) Slater excitations S^'^f, n = 1,2, . . ., with zero eigenvalues. These excitations 
account for the fact that G(3, L + 6) = G(3, L) + 1 [at (3, L + 6) there is one new state, 
/l+6,00) plus Slater excitations of all the (3, L) states] and complete the relative spectrum. 

As an illustration, all the pure relative (3, L) eigenstates for up to L = 17 are enumer- 
ated below: 
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L 


E 


L 


E 


L 







3 


6 


33 
32 ' 


12 


2049 Q Q 
2048 ' ' 


1 




7 


63 
64 


13 


4095 Q 
4UyD ' 


2 


3 
2 


8 


129 
128' 


14 


8193 
8192' "' ^ 


3 


3 
4 


9 


255 n 
256' 


15 


16383 n n 
16384' ^' ^ 


4 


9 
8 


10 


513 
512' " 


16 


32769 Q Q 
32768' ' 


5 


15 
16 


11 


1023 n 
1024' 


17 


65535 n n 
65536' ^' ^ 



Adding a tower of center-of-mass excitations to each of these, one gets the complete 
3-body (LLL) spectrum. 

In the spirit of Eq. (|3l|) , one could also consider "generalized Slater determinants" , of 
the form 



T, 



2ni,2n2,2n3 



s 



= \2n2/~ 



Z2 - 



(36) 



obviously, T2ni,2n2,2n3 times any symmetric function is a bosonic state annihilated by V . 
However, this does not yield new states. E.g., 



422 



35'|/2oo ; 



642 



243 2 r 1 ..4 



(37) 



5 A^-body problem 

As already stated, the diagonalization of the interaction Hamiltonian can be performed 
separately in each (A^, L) sector. An improvement can be made by excluding the center-of- 
mass excitations a priori, i.e., diagonalizing in the pure relative basis — thereby reducing 
the dimension of the subspace involved from G{N,L) to G{N,L). Since the maximum 
possible number of nonzero addends in a partition of L is L itself, G{N, L) ceases to grow 
with N at N = L. Hence, for any A^ > L also G{N,L) = G{L,L) = g{L). Now, by 
definition dZTj), 0(1) = G{L,L) - G{L,L - 1). In each of the G{L,L - 1) unordered 
partitions of L — 1 into L nonnegative addends, at least one addend is equal to zero. 
Replacing that zero with one, we see that G{L, L — 1) is also the number of unordered 
partitions of L into L nonnegative addends of which at least one is equal to one. Hence, 
G{L) is the number of unordered partitions of L into L nonnegative addends of which 
none is equal to one. E.g., ^(6) = 4, which corresponds to the four partitions (zeroes 
dropped for brevity): (6), (4,2), (3,3), (2,2,2). The generating function is 

00 1 00 

n T—I = 1 + E GiL)x^ ■ (38) 

k=l ^ L=l 
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L 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 







1 


1 


2 


2 


4 


4 


7 


8 


12 


14 


21 



Thus, all the pure relative eigenstates for any at a given L can be found by diag- 
onalizing no more than a Q{L) x Q{L) matrix. The basis is formed by the set of states 
{fh...iN '■ J2iLih = L, li > li+i > 0, Zj 7^ 1}. Indeed, any state with any number of 
equal to 1 is a linear combination of the basis states: Since zi - 



-{fl + • • • + fh...h^i 

+ fh...li.i,0,li+l + l-lN + • • • + fh...li-l,0,li+l-lN + l 



(39) 



and if some Ik = 0, then the addend involving + 1 is equal to the LHS, due to symmetry. 
A single k = 1 is thus eliminated, and by repeating, one eliminates them all. 

One pure relative eigenstate for any L and > L is already known, Eq. ()3()|1 . For 
L = 2 and L = 3, as evidenced by the table above, there are no more. For L = 4, apart 
from the state (jSO)), which in this case is 



f(i) 
Jna 



N 



with the eigenvalue 

there is the state 

whose eigenvalue is 

as found in Ref. [T^ . 

For L = 5, the eigenstate 



(1) N{N-?,) 



E 



(2) 
NA 



~ 2 ^ 

N{N -I) 



z,f 



In 



(40) 

(41) 
(42) 

(43) 



N 

^ ^ ^i^jZ^ZiZ^ , 



E 



(1) 

N,5 



N{2N - 7) 



is complemented by a new eigenstate 



N 



■ Z-i -\- Zj 



~Zkf 



(44) 



(45) 



(46) 
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with 



E 



(2) 
TV, 5 



N{N -I) 



15,, 3 
— + - . 
16 4 



(47) 



All the A^-body eigenstates with L < 5 have thus been found exactly. Equations 
and (jTfj) seem to suggest a systematic pattern — which, however, does not exist. 
For L = 6, the basis states are 



hi 
hi 



N 



N 



E~4~2 



N 



i<j 
N 



E~2~2~2 
Z ^ Zi. 



defining a matrix Vij such that 



vhi = J2 v^A 



■J ' 



(48) 
(49) 
(50) 
(51) 

(52) 



one obtains 



( 16Af^-47Af+25 
32 

N+9 
32 

Af-7 
64 





15 
32 

16Af^-60Af+39 



32 

_33 
64 

3Ar+6 



_3 

8 

8Af^-32Af+13 



16 

3 



3 
8 





4Af^-16Af+9 



(53) 



One eigenvalue is 



E 



(1) _ N{N - 4) 



(54) 



which is nothing but Eq. ()30p with L = 6; the other three eigenvalues are the roots of the 
equation 

512^2 + {-768N^ + 2736A^ - 2016)^^ 

+ (384A^^ - 2736A^^ + 6850N^ - 7116iV + 2664)E 

+ (-64A^^ + 684A^^ - 2921A^^ + 6378iV^ - 7527iV^ + 4554A^ - 1080) = . (55) 

For N < L, there are fewer than Q{L) states, but the general scheme of calculation 
remains the same. At L = 6, the rational eigenvalue (j54j) does not exist for < 6, but 
the irrational ones do for all > 4. The table below lists all the pure relative eigenvalues 
for 4 < A^ < 6 and L < 6. 
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L 


= 4 


A^ = 5 


A^ = 6 


2 


4 


15 
2 


12 


3 


3 


25 
4 


21 
2 


4 


9 13 
^' 4 


c; 51 
0, g 


q 21 


5 


3 


15 97 
4 ' 16 


15 81 
2 ' 8 


6 


1.38187, 2.14669, 3.03395 


3.81328, 4.85852, 6.04695 


6, 7.26393, 8.54438, 10.06670 



There appear a few more rational eigenvalues beyond L = Q. In the 4-body problem, 
we already know about the eigenvalue 3 for any odd L, Eq. (jH^ : besides, one of the (4,7) 
eigenvalues is while at (4,12), there duly appears the Slater determinant squared. 



Si = {Zi - Z2f{Zi - 23)^(^1 - ^4)^(^2 - Z3f{z2 - Zif{z3 - Z^f , 

with eigenvalue 0, Eq. Most of the values, however, appear to be irrational. 



6 A model of interacting fermions 

Consider now a Fermi gas with the interaction 



(56) 



which is well-known to have a trivial ground state — the Laughlin FQHE wave functions 
(with the vanishing eigenvalue) [TS|. 

Let us again consider the model projected on the "LLL": it transforms into (again, 
u = l and V/{2n) = 1) 



Zj Zrj 



i<j 



OZi 



Zi,Zj^{zi+Zj)/2^ 



(57) 



Necessarily, one has /(zi, . . . , zjy) = S^hi^Zi^ . . . , zn), where h{zi,...,ZN) is a regular 
A^-body bosonic wave function. In turn, (jHTj) acting on h{zi, . . . , Zn) takes the form 



E 



[_ Zl Zj 



Zi,Zj~*{zi+Zj)/2 



Sn 



{hi 



Therefore, one is back to a situation quite reminiscent of the bosonic problem addressed 
above, and the same machinery applies, but with the new Fermi LLL-interaction term 
given by Eq. (jKH|) . 
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One almost obvious A^-body eigenstate is Sn = Ui^jizi — Zj) [i.e., h{zi, . . . , zn) = !]• 
Indeed, the interaction term (j57|) maps a fermionic state with a given angular momentum 
on a fermionic state with the same angular momentum without introducing any singu- 
larity. It follows that, in the lowest possible angular momentum sector L = N{N — l)/2, 
the single state Sn is necessarily mapped onto itself. This is expressed in the identity 



Sn 



E 

i<j 



Zi,Zj^{zi+Zj)/2 



E 



N 



or 



where 



and 



^ J-J [jZi - Zk) + {Zj - Zk)f ^ ^^^^ 



i<j k^i,j 



[Zi - Zk){Zj - Zk] 



E 



D 



N 



22(N-2) 



D 



N + 1 



N 



'2N + V 
, N , 



is related to the Catalan integer sequence Cf 



Dn = CAN - jH""-'-' 

j=0 



1 (2f 



(59) 

(60) 

(61) 
(62) 

(63) 



The identity can be shown to hold in general by relating Dn-i to the double contour 
integral in the complex plane 



'AT 



dw dz 1 {z + w — 2zky 

(27ri)2 {z - wy {z - Zk){w - Zk) 



such that 



'AT 



2Dn-i + 22(^-1) ATJi 



(64) 



(65) 



and computing In by expanding the contours of integration to infinity [TTj. 
The eigenvalue grows slightly slower than linearly with A^: 





2 


3 


4 


5 


6 


7 


8 


9 


10 


En 


1 


9 
4 


29 
8 


325 
64 


843 
128 


4165 
512 


9949 
1024 


185517 
16384 


424415 
32768 



The asymptotics is En 2N — Ay 

Finally, as examples, here are a few other eigenstates. For all A^: 



(1) 

N 



N 

Sn E 

1=1 



(66) 
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with eigenvalue 



E 



(1) 

N 



En + 



dw dz ^ {z + w — 2zkf' 



22{iv-i) ^ ^2 J (2vri)2 - Zk){w - z^) 



En + 



2N 



227V-2 \ iV - 1 



and 



AT 



f (2) _ Q 

In — 2^ 



zf 



(67) 



(68) 



i=l 



with eigenvalue 



E 



(2) 




dw dz . , -iX- (z + 1(7 — 2zfc)^ 
(z - Zfc)(w; - 2;fc) 



(69) 



Both ()67|) and ()69|) have been obtained by integration in the complex plane, following the 
procedure used in [T7j . 

For = 4, L = 10, one eigenstate is found to be a linear combination of hi = 



i<j\ 



and of the Pfaffian state /i2 = {zi — ,23) (zi — Z4){z2 — -23) (-22 — -24) 



Z2){zi - Z4){z2 - -23) (-23 - z^) + {zi - Z2){zi - Z3){z2 - Zij^Zs - 2:4). Both hi and /i2 can be 
expanded in monomials of Zi as explained in the bosonic case. Denoting fi = S^hi and 
/2 = S4^h2 and introducing a matrix l^^- such that 



Vf^ = E y^^f: 



3 ' 



one obtains 



IK 



/ 43 63 

16 32 

JL 227 

\ 64 128 . 



The eigenvalues are 



(1,2) 
4,10 



571 ±9V393 
256 



(70) 



(71) 



(72) 



Clearly, the machinery developed above in the bosonic case can be thoroughly used for 
the fermionic model considered here. It follows that the A^-body problem is solvable in any 
given angular momentum sector. However, as in the bosonic case, complicated irrational 
coefficients and eigenvalues are expected, as illustrated in a particular case (j72j) . 



7 Conclusion 

Considering the quantum-mechanical model of bosons with a delta-function coupling pro- 
jected on the lowest Landau level, we have completely solved the three-body problem. 
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identified some analytic eigenstates for TV > 4 which belong to two hierarchies (the L — 5 
one is new), and worked out an algorithm through which all other eigenstates can be 
obtained by means of diagonalizing finite matrices (i.e., they are solutions of algebraic 
equations of finite power) . An exact analytic solution of the A'^-body problem is evidently 
out of reach, but a numerical solution to any precision is quite straightforward. A model 
of fermions whose ground state is known to be the Laughlin FQHE wave function, has 
been analyzed along the same lines. We have shown that the Slater determinants, for any 
number of particles, are eigenstates with rational eigenvalues related to Catalan numbers, 
and identified a few excited states. Here too, all the levels are solutions of algebraic 
equations. 
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